Variational formula for experimental determination of high-order correlations of 

current fluctuations in driven systems 
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For Brownian motion of a single particle subject to a tilted periodic potential on a ring, we pro- 
pose a formula for experimentally determining the cumulant generating function of time-averaged 
current without measurements of current fluctuations. We first derive this formula phenomenolog- 
ically on the basis of two key relations: a fluctuation relation associated with Onsager's principle 
of the least energy dissipation in a sufficiently local region and an additivity relation by which 
spatially inhomogeneous fluctuations can be properly considered. We then derive the formula with- 
out any phenomenological assumptions. We also demonstrate its practical advantage by numerical 
experiments. 
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Introduction: Transportation properties in linear re- 
sponse regime out of equilibrium are determined by the 
second-order correlation of time-averaged current fluctu- 
ation in equilibrium [l[. As an extension of this well- 
established formula, a nonlinear response formula was 
derived [2J, where stronger non-linearity arises when 
the third and higher order correlations are more rel- 
evant. Thus, toward a systematic understanding of 
strongly nonlinear transportation, a universal law con- 
cerning high-order correlations must be considered. 

In general, let / be a time-averaged current during a 
time interval r. Its statistical property is characterized 
by a cumulant generating function G(h), which is defined 
by 



rG(h) _ 



(1) 



for sufficiently large r. In recent years, G(h) has been 
studied theoretically for several models [3l4l0l| , partly mo- 
tivated by its symmetry property called fluctuation the- 
orem [ll|. However, it is difficult to experimentally mea- 
sure G(h) when the distribution of I deviates far from the 
Gaussian. Since physically important quantities should 
be obtained experimentally, there should be a formula for 
making G(h) measurable. 

Here, let us recall Einstein's theory for macroscopic 
fluctuation of thermodynamic variables in equilibrium. 
As an example, the cumulant generating function of 
space-averaged magnetization to in magnetic materials 
without a magnetic field is equivalent to the Gibbs free 
energy density as a function of the magnetic field, and the 
thermodynamic functions can be determined by measure- 
ments of heat capacity and susceptibility for any mag- 
netic field and temperature. This implies that the cumu- 
lant generating function of to is obtained without mea- 
surements of fluctuations. A universal idea behind Ein- 
stein's theory is that when a state is determined by an 
extreme condition of a function, fluctuation properties at 
more microscopic scales are described by the variational 
function. The relation between the least action principle 



in classical mechanics and the path-integral formulation 
of quantum mechanics follows the same idea. We thus 
conjecture that the cumulant generating function G(h) 
of / might be identified with a variational function as- 
sociated with a variational principle that determines a 
steady state. 

In this Letter, we study Brownian motion of a single 
particle subject to a tilted periodic potential on a ring. 
Starting from the analysis of a simple example based on 
the least energy dissipation principle we propose a 
phenomenological formula for experimentally determin- 
ing the cumulant generating function of time-averaged 
current. Wc then derive the formula theoretically and 
measure G{h) in numerical experiments. 

Basic idea: We investigate an electric circuit where 
three resistances, Rq, Ri and R2, are connected in se- 
ries under voltage Vq- In particular, in the limiting case 
that Vq — > 00 and Rq — > 00 with / = Vq/Rq fixed, the 
system is assumed to be in the constant-current environ- 
ment. Furthermore, we can impose the constraint that 
the voltages over the resistance R\ and R2 are fixed as V\ 
and Vii respectively. The total energy dissipation ratio 
in this system is given by 



HVi,v 2 \i) 



Ri 



V 2 

R2 



2(Vi + V 2 )I, 



(2) 



where we ignore a contribution independent of V\ and V 2 . 
When we remove the voltage constraint, the values of V\ 
and V 2 are determined by the condition that the current 
is equal to /. Here, these values minimize $(Vi, V2I/). 
This is an example of the least energy dissipation princi- 
ple. 

We next investigate another electric circuit, where two 
resistances, R\ and R2, are connected in series and the 
resistance Rq is connected in a parallel way under the 
total current Iq. In the limiting case that Iq — > 00 and 
Rq — > with V — Rolo fixed, the system is assumed to 
be in the constant-voltage environment. We can impose 
the constraint that the current passing the resistances R\ 



2 



and i?2 is fixed as /. The total energy dissipation ratio 
in this system is given by 



tf(J|V) = (i?i +R 2 )I 2 - 2VI, 



(3) 



where we ignore a contribution independent of /. With- 
out the current constraint, the value of / is determined 
by the condition that the voltage is equal to V. Here, 
this value minimizes ^(/l^). This is another example of 
the least energy dissipation principle. 

Now, a variational principle is expected to be related 
to the description of fluctuation properties as discussed 
in Introduction. For example, the probability of time- 
averaged current fluctuations in the constant-voltage en- 
vironment is given as 



Prob(7|y) ~ e 



,[*(/|V)+co] 



(4) 



where Cq is the normalization constant and the coefficient 
AT respects the fluctuation-dissipation relation and the 
Boltzmann constant is set to unity. Then, the cumulant 
generating function G(h\V) is determined by 

G(h\V) = -±mm[*(I\V)-4ThI + co\. (5) 

The condition G(h = 0\V) = gives c = -%(I S \V), 
where I s represents the steady current. Although ([5]) is 
a simple mathematical expression, we here attempt to 
rewrite it so as to be measured experimentally. The idea 
is to consider a modified system in which an extra voltage 
v in addition to V is applied. We define = vRi/{R\ + 
R2) which represents the allocation of the extra voltage 
v in the circuit. We also denote by I£ the steady current 
for the modified system. By using these quantities, we 
can express ^(I\V) — AT hi + Co as 

2 

*(«!,, V2*\I°) + W - ^ - 2I ( 2Th - V )- ( 6 ) 
i=l 

By choosing v — 2Th and minimizing in terms of /, 
we obtain G(h\V) = -$(<ui*, v 2 *\I%) / {AT) . The varia- 
tional principle in the first example leads us to rewrite it 
as an additive form 



1 2 
G{h\V) = min V 



11 2T v v- 

Rt 



(7) 



This additivity relation is the basic formula which we ex- 
tend so as to apply more general cases. Since the voltages 
can be controlled experimentally, the right-hand side of 
(UJ) is obtained experimentally. Although the formula is 
derived for the simple system, the result is quite sugges- 
tive. The assumption for the formula is that the least en- 
ergy dissipation principle and the fluctuation-dissipation 
relation hold locally in space. Although this assumption 
is not always valid, many non-equilibrium systems are 
included in this class. Below we apply the relation ([7]) to 
the Brownian motion out of equilibrium. 



Main Result: The system we study consists of a single 
Brownian particle on a ring with a size of L. Its motion 
is described by a Langevin equation 



7± = / 



dU 
dx 



(8) 



where x = dx/dt and £ is the noise satisfying (£(i)£(t')) — 
2T^fS{t — t'). 7 represents a friction constant, / a 
uniform driving force, and U{x) a periodic potential. 
The external force acting on the particle is written as 
F{x) = f — dU{x)/dx. It should be noted that the 
Langevin equation ([8]) corresponds to experimental sys- 
tems that have been studied with the aim of tes ting 
new ideas in non-equilibrium statistical mechanics 12- 
14[. We are particularly concerned with the cumulant 



generating function G{h) defined in ([1} with / replaced 
by 



1 



dx 
dt—. 

dt 



(9) 



Throughout this Letter, the boldface font is used to in- 
dicate a function of x. 

Here, in order to find a formula similar to ([7]). we apply 
an extra force u{x) in addition to F{x) with an interpre- 
tation that F{x)Ax and u{x)Ax correspond to the volt- 
age and the excess voltage over the interval [x,x + Ax], 
respectively. We consider an empirical probability den- 
sity p{x) and its current j{x) obtained by measure- 
ment during a finite time interval r. That is, p{x) = 
f Q dtS{x — x{t))/r and j(x) = j Q dtS{x — x(t)) o x/t, 
where o is the Stratonovich rule of the multiplication. In 
the limit r —t 00, j(x) — > J" and p(x) — > p™(x), where 
the steady probability density and the spatially 

homogeneous current J" for this modified system with 
the extra force u are determined by 



J? = ±p?(F + u)--d x pl 

7 7 



(10) 



This shows that the resistance in the interval [x, x + Ax| 
corresponds to r y(Ax)/p™(x). Finally, since the current / 
corresponds to dxj(x)/L = v/L, h in the formula for 
the time averaged current is replaced by hL for the case 
of cumulant generating function G(h) for v. That is, the 
constraint condition v\ + v 2 = 2Th becomes 



/ dxu(x) = 2ThL. 
Jo 



(11) 



From these correspondences, we heuristically derive 

G(h) = ~ min g{u), (12) 

4J u;f£ dxu(x)=2ThL 



with 



Q(u) = / dx 



^u(x) 2 



2J>(x) 



(13) 
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That is, G{h) is proportional to the minimum energy 
dissipation of the extra force u(x) when the system is 
assumed to be in constant-current environment of the 
current J s u with J Q L dxu(x) = 2ThL. The formula (12) 
with (JTSJ) is the main claim of this Letter. 

We determine the optimal force u op t(x) that satisfies 
G(u op t) = ini u Q(u) with the condition (jlll) . By calcu- 
lating 5Q = Q(u + Su) — Q{u) explicitly, we find that 
SQ = leads to 



K=-{ F 

7 



^opt 



^opt 



— Or, 



'a;^opt i 



(14) 



where K is a constant. Differentiating the both-sides 
with respect to x, we obtain u op t(x) under the periodic 
boundary condition. The substitution of the obtained 
result u op t into (|14l) leads to K. This constant K is 
significant. Indeed, by multiplying p s opt to (|14p and u op t 
to (|10|) and by considering the difference between them, 
we obtain K — 2TG(h). Furthermore, for the non- linear 
differential equation (TUT) with K = 2TG(h), we perform 
a Cole-Hopf transformation defined by u opt (x) — 2T(h + 
d\ogip(x)/dx), where ip( x ) > 0- We then derive a linear 
eigenvalue equation Chip = G(h)ip, where the operator 
Ch is defined by 



C h - = -(d x + h)-+-{d x + h) 2 

7 7 



(15) 



From the positivity of if), G{h) turns out to be equal to 
the maximum eigenvalue A max of the operator Ch and 
w op t is related to the corresponding eigenfunction. Ex- 
amples of G(h) thus determined uniquely are displayed 
in the left-side of Fig. [TJ 

Now, we prove (|12p without employing any phe- 
nomenological assumptions. For a stochastic variable 
y(t) = tv(t) = J Q dsdx(s)/ds, we denote by P(y,t) the 
probability density that y takes the value y at time t. We 
then have exp(rG(/i)) = J dyP(y,r) exp(hy). Here, the 
function Q(y,t) = P(y,t) exp(hy) obeys d t Q = MhQ, 
where 



1 



T 



(16) 



M h - = -~(d y ~ h)(F(y)- ) + ~{dy ~ hf ■ 

7 7 
We denote the maximum eigenvalue of the operator Ait 
by Mmax- Since Q(y,t) ~ e' lmax ' for sufficiently large 
t, G(h) is equal to fi max . (See a related result in Ref. 



[9j.) By comparing (|T5| and (|T6| . we find L h = Mh- 
Thus, A m ax = Atmax. Since the eigenvalue equation 
Chip = G(h)ip is equivalent to the variational form (fT2|) 
with (fl3l) . we end the proof. This method of derivation 
of the formula can also be applied to Langevin equations 
in any dimensions. See Ref. [15j for details. 

The formula (fT2]) is closely related to the so-called ad- 
ditivity principle [3J47J. We mention the relation explic- 
itly. Let us consider the probability density of J = v/L, 
which is expected to possess a large deviation property 



Prob( J = J) ~ c 



-rr(J) 



(17) 



Through mapping from ([8]) to fluctuating hydrodynam- 
ics [lj| , we can apply the additivity principle to the con- 
tinuum description of the fluctuating density field. The 
result is 

r-L 



r( J) = min 
p 



dx 



(J-f s (x;p)Y 



4Tp(x)h ' (18) 
u(x))/^(—Td x p/^ for any u(x). 



where j™(x;p) = p{F(x 
Since G(h) is connected to T(J) as G(h) — maxj[hLJ 
r(J)], we obtain 

r-L 



G(h) 



1 . 
mm 

4T J, P 



where 



p(x) 



dx(f)(x; J, p), 



u{xf - 2u(x)j?(x; P ) + ^-AJ- P) 



(19) 



7 p{%) 

under the constraint condition (|11[) . We make a variable 
transformation from p to it by d x j™(x; p) = 0. Explicitly, 
for a given u, we can determine J" and p™(x) uniquely 
as j™(x;pf) — J". The minimization with respect to 
J in (fT9|) is achieved by J = J", and / Q L dxcf>(x; J, p") 
becomes t/(n). Thus, (fT9|) leads to the formula (fT2|) . 

Experiment: We consider the experimental measure- 
ment of G(/i). In the standard method, G(/i) may be 
determined from the measurements of cumulants. How- 
ever, high-order cumulants are quite difficult to be mea- 
sured experimentally. Here, we claim that our formula 
can provide G(h) experimentally if we can control the 
external force and know the value of the temperature T 
and the friction constant 7. The method is as follows. 

We measure a trajectory (x(i))o<t< T for the system 
with an extra force w(x). From this data, we define 

~w{x{t)) 2 



Gr{w) = 



1 



dt 



7 



2x(t) o w(x(t)) 



(20) 



Note that G T (w) — > G(w) in the limit r — > 00. Thus, if 
we know the optimal force it op t , the approximate value of 
G(h) is obtained from Gt(w) by setting w ~ u op t- Here, 
we have an identity 

G{u opt ) = -AThJ™L - [ dx^u opt (u opt - 2w) (21) 
Jo 7 

for any w, which can be derived from the calculation 
of G(u opt ) -G{w) with the aid of (fTTJ|) and ([H]) [l5|. 
We express u op t(:c) = 2Th + J2n=-N where 
4>n{x) = cos(27rnx/L) for n > 0, <p n (x) = sm(2Tmx/L) 
for n < 0, ao = 0, and N is the truncation number 
of the approximation of the force u op t(x). Since (|21j) 
holds for any w, we determine {a n )^ = _ N by preparing 
2N + 1 forces w^(x), —N < i < N. Concretely, by us- 
ing the trajectory x^ l \t) for each w^ l \x) = 5T4>i(x)/L, 
the right-hand side of (f2~Tj) is approximated as 



H T (w 



T 



dtU {€) {t), 



(22) 



4 



where 

H^(t) = -AThx^(t) + 2^ n M^(t))w^{x^(t))a n 
E B ,m ( x{i) (*))&» (* W (t))a n a m . (23) 

Since Hoofa^') is independent of ujW, we consider the 
following equation for (a n )n=-N' 

H T (w {l >) = H T (w<- t+ ^), (24) 

where — N < i < N. By solving (|24p. we obtain an 
approximation of M op t(x). The approximation becomes 
more accurate for larger N and larger r. In the right-side 
of Fig. [TJ we display an example of the measurement of 
G(h) in numerical experiments. 
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FIG. 1: (color online). G(h) for the Langevin equation © 
with U(x) — Uo cos(2irx/ L). Quantities are converted to di- 
mensionless forms by setting 7 = T = L = 1. We fix / = 1. 
(Left:) Numerical calculation of K = 2TG(h) in ff^ for 
Uo = (red solid line), 3 (green dashed line), and 5 (blue 
dotted line). (Right:) Experimental measurement for the 
system with Uo = 3. We assume to know T = 7 = 1 and 
determined G(h) experimentally from trajectories (x(t))J=o 
according to the method described in the text. N = 5 and 
r = 250 (triangle symbols), 1000 (square symbols), and 4000 
(circle symbols). The error-bars are within the symbols. The 
dashed line is the same as that in the left-side. 

Concluding remarks: In this Letter, we have proposed 
the formula (fT2"|) for a driven Brownian motion described 
by ©. We have demonstrated that G(h) can be deter- 
mined experimentally from T, 7, and trajectories of the 
modified system with an external force. Note that the 
formula yields an expression of the fe-th order cumulant 
Cfe = T k ~ 1 (v k ) c 15], where C\ and C2/2 are equal to the 
average velocity and diffusion constant of the Brownian 
particle, respectively. The result reproduces the known 
expression for the diffusion constant 17j in the simplest 
manner. Furthermore, (|12j) can be extended to more 
general bulk-driven systems when we assume fluctuat- 
ing hydrodynamics. We express G{h) by a path-integral 
expression, and if it is determined by the contribution 
of a saddle configuration, the formula (fT2"j) is valid, as 
for the case of (J5J. However, there are many systems 
for which the condition is not satisfied, as investigated in 
Refs. 4-7]. (See also sections 3 and 4 in [l8[ for directly 



related discussions on the problem in this Letter.) Future 
work will concern the determination of the range of the 
application of the formula. We believe that our theory 
for the derivation provides a sound approach in studying 
this problem. We also wish to present non-trivial predic- 
tions by employing the formulas for interacting particles 
systems. We hope that our formula will be studied the- 
oretically, numerically, and experimentally, so that the 
novel nature of non-equilibrium fluctuation can be un- 
covered. 
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